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Abstract 

Dynamical Phase transitions (DPTs) describe the abrupt change in the dynamical proper- 
ties of open systems when a single control parameter is slightly modified. Recently we found 
that this phenomenon is also present in a simple model of a linear array of metallic nanopar- 
ticles (NPs) in the form of a localized-delocalized DPT. In this work we show how to take 
advantage of DPTs in order to design a new kind of plasmonic sensor which should own some 
unique characteristics. For example, if it were used as plasmon ruler it would not follow the so 
called universal plasmon ruler equation [Nano Letters 2007, 7, 2080-2088], exhibiting instead 
an on-off switching feature. This basically means that a signal should only be observed when 
the control/measured parameter, i.e. a distance in the case of plasmon rulers, has a very precise 
and pre-determined value. Here, we demonstrate their feasibility and unique characteristics, 
showing that they combine high sensitivity with this on-off switching feature in terms of dif- 
ferent distances and local dielectric constants. This property has the potentiality to be used in 
the design of new plasmonic devices such as plasmonic circuits activated only under certain 
environmental conditions. 

Keywords: plasmonics, one dimensional arrays, open systems, localization, plasmon rulers. 



INTRODUCTION 

The study of localized surface plasmon (LSP) resonances over the last couple of decades has 
led to great advances in several areas of science and technology. Probably, its most significant 
application is its use in nanosensors.-"— Most of the strategies used to turn a plasmonic device into 
a sensor are based on the shift of the resonant position of the LSP-~— or in exploiting hot spots to 
increase spectroscopic signals.—"— In this respect, there is a growing interest in using interacting 
plasmonic systems to measure distances at the nanoscale. These devices are known as nanorulers 
and they are finding interesting applications in several fields such as biology.™ Most of the proposed 
plasmonic rulers work by the same principle, the shift of a resonance induced by the variation of 
a distance. As noticed originally by Jain et. al.,^ the frequency shift of the LSP resonance follows 
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an almost universal exponential or quasi-exponential law with respect to the separation between 
NPs. This behavior is known as the universal plasmon ruler equation (UPRE) and was observed 
in simulations and experiments using NPs of different size and shape,-"— core-shell NPs,-l^ and 
even two dimensional arrays of NPs.— There are some deviations of the UPRE due to retardation 
effects or the effects of higher order multipoles, but they can usually be corrected by empirical 
fittings.- Only when NPs are almost touching each other, there are some deviations, nonfitable to 
exponential, of the UPRE.— There is another kind of plasmon rulers that is now being actively 
investigated. They are based on a different principle, the fluorescence resonance energy transfer 
(FRET). Its working principle is the resonant transference of excitation from a fluorescent molecule 
to a nearby NP. The closer to fluorophore to the NP, the greater the photoluminescence quenching. 
This effect has been successfully used as a nanorulers in several experimental examples.—""— 

Collective effects induced by LSP couplings in NPs arrays can give rise to other phenomena 
and new potential applications.—"^ One example is the use of Fano resonances to design plas- 
monic sensors.^ In this case, the shift of very sharp peaks or valleys in the extinction spectrum, 
consequence of collective effects of interacting NPs, have been proposed to measure distances or 
the dielectric constant of the surrounding medium.^ More exotic applications of NPs coupling are 
also possible and many of them are now being actively explored, such as clocking by metamaterials 
or plasmonics circuitry.—"^ 

Although, linear arrays of metallic NPs have been proposed as potential waveguides to transfer 
massive amounts of information at the nanoscale,— nowadays it is clear that high damping factors 
would impose severe limits to that and most of the effort in this field is being put on overcoming 
this limitation. However, one important and very general feature of 1-D systems, which seems 
to have been overlooked in plasmonics, is the extreme sensitivity of their dynamical properties 
to slight changes of their parameters. In this respect, we have previously studied the plasmonics 
energy transfer from a locally excited NP (LE-NP), the 0th NP in the scheme shown in[Tl to the 
interior of a semi-infinite NP chain.— This system presents a form of dynamical phase transition 
(DPT), a phenomenon that is currently attracting interest as source of novel effects on various 
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fields.— — Basically a DPT in this context means that by sightly moving a single parameter, such 
as the shape of the NPs or the distance between them, the system undergoes an abrupt transition 
from transferring all the injected excitation to the interior of the NP array (a delocalized state), to 
keeping all the external excitation in the NPs closer to the excitation point (a localized state). ^- 
The aim of this work is to exploit this basic property into a realistic setting, in order to design a 
new kind of excitation-transfer plasmonic nanosensor (ETPS). 

1 2 ... 00 

Figure 1 : Scheme of the system studied. An external source locally excites the LSP of a NP which 
is coupled to a waveguide where detection takes place. The arrow on the right indicates that the 
array is infinite in this direction. 



THEORY 
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The system studied, depicted m\Tl is modeled using the coupled dipole approximation 
In this model, each /*-NP is described by a dipole Pi induced by the electric field produced by the 
other dipoles, Ei j, and the external source, E^^^^\ We assume a generic ellipsoidal shape for the 
metallic NPs and describe its polarizabilities in a quasi-static approximation,— resulting in: 

where Y is the volume, L is a geometrical factor that depends on the shape of the NP and the di- 
rection of E, Eq is the free space permittivity, and £,n,i is the dielectric constant of the host medium 
around the i-th NP. The dielectric constants e, of the NPs are described by a Drude-Sommerfeld 
model: e = £cx, — , where (Op is the plasmon frequency, 7] is the electronic damping fac- 

tor, and £00 is a material dependent constant that take into account the contribution of the bound 



electrons to the polarizability. The LSP oscillations over each NP can only be transverse (T) or 
longitudinal (L) to the axis of the linear array and we assume, for simplicity, that the ellipsoids 
axes are aligned with respect to these T and L directions. If the excitation wavelength is large 
compared with the separation between NPs, Eij can be taken in the near field approximation, 

_yT.Lp, 

^ij = A — where 7 = 1, 7^ = — 2 and d is the distance between NPs. Taking into account 
these considerations. Pi and E^^^^^ can be arranged as vectors P and E resulting in:— 

P= (I«2-M)"^ME, (2) 



where M is the dynamical matrix and M is a diagonal matrix that rescales the external applied field 
according to local properties: 

R,, = -eor^co^fi, (3) 

Because of the cubic dependence of E on d, we neglect contributions beyond nearest neigh- 
bors, ^^'^^ resulting in a tridiagonal M. For NPs belonging to the chain, (z, j) 7^ 0, M,- = co^^ — irjco 
and Mi J = (o^, and for the LE-NP, Mo,o = (Ojpo - irio(o, Mo,i = (o^^ ^ and Mi^ = 60^^. The LSP 



frequencies co^^. and coupling constants co^. ., result: 

T (olLi 
co^ = — (4) 

^'^ 47t£,ndfj 

For spherical NPs (L = 1/3) in vacuum (£,„ = 1) within a Drude-Sommerfeld model which com- 
pletely neglects contribution from bound electrons = 1), Eqs. ??-?? reduce to the simple form 
used in several references.— "^i^^*^ This shows the meaning of //, which essentially collects the 
deviations from the ideal model of spherical NP in vacuum with £cx, = 1 , not only for Ri^i but also 
for couplings between NPs. Note that these deviations reduce to a constant for £00 ~ £m,i and/or 
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small excitation frequencies (strictly for [(0^ — /WT],] ^ w^/fCoo — £.m,i\)- The form of Eqs. ??-?? 
not only extends our previous expressions^^ to non-spherical NPs but also includes explicitly the 
interdependencies of G)^. ., CO^p., and Ri i with the parameters of the system, i.e. size, shape, and 
material of NPs, distances between NPs, and dielectric constants. It should also be mentioned that 
although this model can be further improved, the different corrections will only add quantitative 
corrections as long as the essential physics of the problem remains, i.e. a semi-infinite chain with 
nearest neighbors interactions.— "^2*^ We will analyze this point in more detail in the last section. 

In the present work we will consider the case of excitation injection to a semi-infinte homo- 
geneous linear array of NPs, where solely the LE-NP can be different from the rest and the only 
nonidentical separation is that between the LE-NP and the first NP of the chain. Therefore, there 
will be two LSP resonances, that of the LE-NP co^^^^ and that of the chain's NPs O)^^. Similarly, the 
couplings between the chain's NPs will be all the same ft)^. The couplings between the LE-NP and 
the first NP of the chain are (0^^ ^ and (O^i o' They are not necessarily equal (see Eq. ??), since they 
describe quite different situations, i.e. how P\ affects Pq (co^^ , ) and how Pq affects Pi {co^^ ,,). 

A.U,1 A.I,U 

The values of the resonance frequency co^^ and the couplings co^.. can be tuned in different 
ways. However, one must be careful because they are always interrelated. For example, choosing 
a different material for the NPs will change both ft)„ and Cc^ and therefore co^ and m . Changing 
the shape of the NPs will modify L and consequently O)^ and ft) , except for e^^ = e„, in which 
case it will only modify (0^^. The volume and distances between NPs will change (O^... If the 
material, size and shape of NPs are fixed, the free parameters that can be used for sensing purposes 
are the dielectric constant of the medium £,„ that alter both co^^. and co^.., and the distances between 
NPs that change couplings. 

The values of co^ and co^^ define the passband which is the interval of frequencies at which an 

excitation can propagate with only a relatively small decay, given by rj , through the waveguide 

formed by the linear array of NPs. All excitations with frequencies CO outside the passband will 

decay along the chain exponentially and very fast. In the weak damping limit (WDL) T] — t- 0, this 

passband is given by a simple expression, ^ ^ 2. This fact shows one of the roles of the 

% 
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chain: to determine which frequencies can be propagated. The other role is to perturb the local 
density of plasmonic states (LDPS) of the LE-NP by pushing the resonance from its uncoupled 
value at CO^py to the nearest edge of the passband.-^^*^'^ This process occurs very abruptly when 
the resonance is about to cross the passband edge that is when the dynamics of the system changes 
completely .-^^^ 

This phenomenon can be completely understood in terms of the divergences, or poles, of Eq. 
?? in the WDL.— Basically there are three different dynamical regimes of interest for excitation 
transfer: 1) delocalized or resonant state regime, where the maximum of LDPS of the LE-NP falls 
within the passband and most of the excitation can be transferred to the NP array; 2) localized 
state regime where there is a very sharp peak in the LDPS outside the passband and most of the 
excitation remains close to the excitation point; and 3) a transition regime called virtual state where 
the LDPS presents a non-Lorentzian asymmetric peak just at the passband edge and most of the 
excitation can be transferred to the waveguide. The transitions between the different dynamical 
regimes can be obtained analytically in the WDL and assuming frequency independent couplings, 
where:— 

2/3+4\/2_/32_i 

« = ' 

gives the resonant-virtual transition and: 

a = 2±(i^. (7, 

gives the vz>?Ma/-Zocafee J transition, where a — (co^j ^(O^^^/g)^, /3 = {(o^^^/ (o^^), and V = (co^/co^p). 
The constant a is the relative coupling between LE-NP and first NP of the chain in units of the 
coupling in the chain, j8 yields the difference in resonance frequency of the LE-NP with respect to 
that of a chain's NP, and V measures the strength of the coupling between chain's NPs in units of 
their own resonance frequency. These formulas, although strictly valid only in the WDL, give an 
excellent estimation for finites T]s.— 

The excitation transferred from the LE-NP to the m-th NP of the chain can be calculated exactly 
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(in the quasi-static limit and within the nearest-neighbors approximation) by using the formula:— 

= ""-^f"'^' c->-ii^^^^ (8) 

0)2 ^ [cq2 - (o|po] - an((o)) 

where Eq is the electric field on the LE-NP and the renormalized eigen-frequency Wgp is given by 
G)^ = (0^ —irjOJ. The decay length / and the wavenumber q depend on the self energy n as: 



l±iq = ln{(ol/n), (9) 



where IT is: 



n((0) = i[(02_^2j_ (^Q) 

sgn((o2 _ coiy^^[co^-alf-4col 

In all the plots of the next section with use \Pm\'^ as an indicator of the excitation transferred 
to the waveguide, with w^/w^p = 0.45, rj / co^^ = 0.03, and using m = 8. We have chosen m = 
8, because around this value the behavior of \Pm\^ vs m becomes essentially a small decaying 
exponential for all CO within the passband or zero for excitation frequencies outside it. Thus, 
choosing a bigger value of m would only add a multiplicative factor to the results, i.e. it would not 
change qualitatively the figures shown. 



RESULTS 

In this work, our purpose is to show how to use the sudden change in the dynamical properties of 
systems tuned around their DPT to device ETPSs. For this purpose, we used the system depicted 
in[T]as an example. |2] shows the maximum excitation transferred to a NP of the waveguide, |/'8p, 
enabled by a variation of the excitation frequency CO, at each system configuration. Superposed are 
the critical values that separates dynamically distinct behaviors calculated in the WDL (Eqs. ?? 
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Figure 2: Maximum value, allowed by a variation of the excitation frequency, of the excitation 
transferred to a NP of the waveguide as function of the system's parameters, gray lines indicate 
DPTs: the continuous line corresponds to the resonant-virtual DPT and the dashed lines show the 
virtual-localized transitions. 

and ??). As it has been previously reported,— excitation transfer is controlled by the DPT where 
virtual states are transformed into localized states (gray dashed lines). 

Analyzing the phase diagram of [2lit is clear that the strong dependence of excitation transfer 
to the system's parameters can be used for different forms of sensors. As a first application of 
this, lets us analyze the potentiality of this plasmonic device as a plasmon ruler. In [2] a vertical 
displacement corresponds to a variation of the distance between the LE-NP and the chain. The 
vertical blue arrow indicates an example of a region in the parameters space that can be used for 
such a purpose. For each point along this arrow the system depicts a different excitation transfer 
spectrum. Three representative points along this line, at do-i/d = 0.8, 0.9, and 1.0, are used in 
|3]-a) to illustrate how the spectrum of excitation transfer changes when the system undergoes a 
DPT. Note the abrupt change in the spectra induced by a slight change of doi/d: the excitation 
transfer suddenly increases at (O^/ co^^ = 0. 1 when Jqi / varies from 1 .0 to 0.9, and then decreases 
equally abruptly when doi/d change from 0.9 to 0.8. This behavior is a consequence of the virtual- 
localized DPT, indicated in [21 by the curved gray dashed lines and predicted by Eqs. ?? and ??. 
One of the interesting aspects of this new form of plasmonic ruler is that the variation of a distance 
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Figure 3: a) Excitation transferred to the 8th NP of the waveguide, iPgp in arbitrary units, as a 
function of O) for three different values of do-i/d. Each spectrum corresponds to a point along the 
vertical blue arrow shown in [21 The value of l/'sl^^t ^ p^^^) in El corresponds here to the maximum 
value of P in a given curve. The two vertical grey lines mark the passband edges of the waveg- 
uide in the WDL. b) The same but moving the SFO for fixed (do-i/d) along the horizontal green 
arrow shown in [21 

not merely modifies the position of a peak but dramatically changes the whole spectrum, and quite 
remarkably, before and after the DPT, only negligible signals should be observed independently 
of the excitation frequency. In this sense, our excitation transfer nanoruler could be used as a 
robust nanometric standard. This is so, because, in order to observe a difference in the excitation 
transferred, it is not required a high accuracy in the excitation frequency nor a monochromatic 
light but keeping the maximum of the excitation frequency close to the lower band edge frequency, 
in our example at 0)^/o)^^ =0.1, is a sufficient condition. Furthermore, if the polarization has 
two components one transversal, T, and one longitudinal L, this should not interfere in a sensor 
tuned to work with only one of them. This is so, because the passband of T modes are half of L 
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modes, Eq. ??, which imply that regions in the spectrum where DPTs can be observed, should be 
well separated. Moreover, according to Eqs. ?? and ??, DPTs should occur for different values 
of system's parameters. Another important aspect is that the detection could be performed in a 
region encompassing several NPs. This is so, because, except for the first NPs of the chain, the 
excitation transferred to the chain is an exponential function of NP position m, see Eq. ??. Thus, 
the main feature of this kind of nanorulers is its robust on and off feature, i.e. signals should only 
be detected when the LE-NP is at a very precise and predetermined distance. 




Figure 4: a) Relative increment of the excitation transferred, {\P\-^ — \P\j„itiai)/\P\hiitiai' ^ func- 
tion of the relative expansion of the whole system, Ad/dinUiah for w^/ co^^ = 0. 1. The initial con- 
ditions are do^/d = 0.97 and SFO=-1.58. b) The same but changing the local dielectric constant 
around the LE-NP, £^,0- The initial conditions are do^i/d — 0.83 and SFO=-1.33. 

ETPS could also be used for sensing other properties, not only as plasmon nanorulers. As 
shown in [2l the system is also very sensitive to the value of the square frequency offset (SFO = 
(^s^po ^ ^sp) / consequently to CO^^^, CO^^, and/or co^. Increasing the numerator of the SFO 
or decreasing its denominator, will both move the SFO away from the center of|2]as indicated by 
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the horizontal green arrow. [31-b) shows the excitation transfer spectrum for three values of SFO 
along this green arrow. Similarly to the nanoruler case, the spectrum experiences a significant 
change when moving the SFO along a DPT, from SFO = —1.5 to —1.8. In order to show how to 
make profit of this feature we will consider two examples: one in which contraction or expansion 
of the whole system is measured, and one in which local dielectric constant around the LE-NP is 
measured. |4]-a shows, for a fixed excitation frequency, £0^/ co^^ = 0. 1, the relative increment of the 
excitation transferred, — \P\j„itiai)/\P\fnitiai' a function of the relative expansion/contraction 
of the whole system, which turns our device into an optical strain monitor^ or even into a temper- 
ature sensor, depending on the expansion coefficient of the host material. Note that a relative 
expansion/contraction of the whole system corresponds to a change in co^ but not in do-i/d, CO^p^y 
nor in co^^. Thus it is equivalent to an horizontal displacement in [21 IH-b shows the variation of the 
relative increment of the excitation transferred as a function of the local dielectric constant around 
the LE-NP, a property potentially useful for molecular sensing purposes for example.^ In this case, 
a variation of £,n,o changes not only (U^p^ but also £0^^^^ and Rq q (Eqs. ??-??). Therefore, it is 
equivalent to a diagonal displacement in [2] and a rescaling of the external applied field. 

RETARDATION EFFECTS AND CONTRIBUTIONS BEYOND 
NEAREST NEIGHBORS 

The key ingredients from which the phenomenology described so far arise are: dominant nearest- 
neighbors interactions and the semi-infinite character of the system. If this conditions are fulfilled, 
the results presented above should be applicable to the system under consideration besides quan- 
titative corrections such as shifts of the spectra or the exact position of DPTs in the parameter's 
phase space. However, it is important to understand how our results are corrected by different 
factors. 

In this section, we will study the effect of interactions beyond nearest-neighbors and retardation 
effects. In this case, Eq. ?? is still valid but all coupling between NPs should be considered, i.e. M 
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is no longer tridiagonal, and these couplings should be determined by the dipole-induced electric 
field beyond the quasi-static limit, i.e.: 



{{kdfidx p)xd+ [3d{d-p)-p] (l-ikd)}, (11) 



where k is the wavenumber in the dielectric, k= (o/v (where v is the speed of light in the medium), 
d is the unit vector in the direction of d (where d is the position of the observation point with 
respect to the position of the dipole). 

As the system consists of a linear array of NPs where the spheroids axes are aligned with 
respect to the direction of the array, transversal {T) and longitudinal (L) excitations do not mix. 
Thus, the coupling terms a)| of Eq. ?? acquire the form: 

<=^'^^fi^ (12) 



where 



%j = -2[l-ikd,,jy'''--J 

%^ = [\-ikd,,j-{kd,,jfy^'''-K (13) 

Eq. ?? has now to be solved numerically which can be done by using standard methods of linear 
algebra.— 

The physical system considered in this section corresponds to a linear array of very flat oblate 
spheroidal NPs of silver (Coo = 3.7, and (Op = 9.3xl0^^rad/s) with the radio of the minor axis equals 
to 4.5nm and the radii of the major axis equal to 46.5nm (shape factor L, = 0.86,). The major to 
minor axis ratio is about 10 which is large enough to ensure a strong quadrupole quenching.— 
Separation between NPs of the chain is 13.5nm, from center to center, and they are aligned with 
the two equal axes perpendicular to the array direction. The external electric field is applied only 
to the first NP where the direction of Eq^'^ is also perpendicular to the array, i.e. J =yJ j. Elec- 
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tronic damping is calculated using the Matthiessen'srule 7] = v/r(l/Zbuik + C'//eyy) where the Fermi 
velocity is Vf = 1.38xl0^m/s, the bulk mean free path is /bulk = 51nm, we use C = 1, and the ef- 
fective mean free path is calculated using Zeff = 4-V /S where V is the volume and S the surface of 
the spheroids. For simplicity we consider a dielectric medium with £,„ = e^. Taking into account 
this, T?/ft)gp and tt)^/w^p give 0.032 and 0.41 (for kd = 0) respectively . We use here a slightly 
smaller coupling than that of the previous section for reasons that will become clear in the context 
of [5l As before, the observation point is fixed at the 8th NP, and we use the value of \P%\^ as an 
indicator of the excitation transferred to the chain. Due to the large interparticle coupling, finite 
size effects result so important, especially near the lowest band edge, that it is necessary to take 
very long size of chains. We use 200 NPs, a value large enough to ensure negligible finite size 
effects. 
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Figure 5: Excitation transferred to the waveguide, iPgP in arbitrary units, as function of (0^/ (O^^ for 
dQ-i/d = 0.9, «x/®sp = 0-41 and T7/fi)^p = 0.032. Different curves show numerical results within 
the near field approximation for: nearest-neighbor interactions, fifth-nearest-neighbor interactions, 
and a fully converged calculation. The system is at the DPT for the curve with nearest-neighbor 
interactions. 

[5] shows the effect of including interaction beyond nearest-neighbors. The first effect of this 
is essentially a red-shift and a broadening of the passband. The second effect is that the DPT 
is also shifted from do-i/d = 0.9 to do-i/d = 0.84 as can be seen in[6l-a. However, it is still 
present and indeed, as the damping term is proportional to the frequency, the peak is increased. 
Note that if (O^ / co^^ = 0.45 were used, the peak would disappear from the spectrum as it would 
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be located at co^ < 0, showing an important effect of higher-order contributions to interactions. In 
[6l-a and b retardation effects on the DPT are compared. Here, retardation causes nothing else but 
an increase of the intensity of the peak, in the region of the spectrum where the DPT is observable. 
This is because, in this region of the spectrum, around 0)^/ 03^^ = 0.02, d is much smaller than k 
(kd ~ 0.05). Of course at higher frequencies, for example at (O ^ cosp, this is not the case (kd ~ 0.4) 
and the behavior of the system is completely beyond the near field approximation. 
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Figure 6: a) Numerical results without retardation effects of the excitation transferred to the waveg- 
uide, l^gp in arbitrary units, as function of aP'/co^^ for different values of do-i/d. b) The same 
but with retardation effects. In both cases all interactions between NPs are considered. 

Finally, |7] shows the value of IP^I^ for different NPs of the chain. Note that even after including 
all contribution to the interaction among NPs and retardation effects, \Pm\^ still decays exponen- 
tially as predicted by the theory, although the decay rate is different mainly because the passband 
edge is considerably shifted. 
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Figure 7: Excitation transferred to a NP of the waveguide, \Pm\^ in arbitrary units, as function of 
its position m. Retardation effects as well as all interactions between NPs are considered. Here 
do-i/d = 0.84 and the excitation frequency is (0^/(0^ = 0.02, which correspond to the peak of 

CONCLUSIONS 

In this work we have studied the potentiality for sensing purposes of the excitation transferred 
from a LE-NP to the interior of a semi-infinite NP chain, when the system is very close to its 
DPT condition. While most plasmonic sensors use the shift of the LSP resonance peaks or the 
photoluminescence quenching of hybrid plasmonic-fluorophores systems, here we introduce a new 
working principle. Basically, the idea is to take advantage of the abrupt change in the plasmonic 
energy transferred when a control parameter is slightly changed around a DPT. We have shown 
that this kind of sensor has the unique characteristic of having an on-ojf switching property and 
a high sensitivity, which opens new possibilities to design plasmonic devices such as plasmonic 
circuits activated only under certain environmental conditions. We have also addressed the effects 
of different corrections to the approximations used and have shown that if the system's parameters 
are chosen properly, the DPT should be observable and therefore useful for the sensing purposes. 
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